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Abstract 

The aim of this paper is to construct a natural Riemann-Lagrange 
differential geometry on 1-jet spaces, in the sense of nonlinear connec- 
tions, generalized Cartan connections, d-torsions, d-curvatures, jet elec- 
tromagnetic fields and jet electromagnetic Yang-Mills energies, starting 
from some given nonlinear evolution ODEs systems modelling biologic 
phenomena like the cancer cell population model or the infection by hu- 
man immunodeficiency virus- type 1 (HIV-1) model. 
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1 Historical geometrical aspects 

The Riemann-Lagrange geometry [5] of the 1-jet space J^(T, M), where T 
is a smooth "multi-time" manifold of dimension p and M is a smooth "spa- 
tial" manifold of dimension n, contains many fruitful ideas for the geometrical 
interpretation of the solutions of a given ODEs or PDEs system [7]. In this 
direction, authors like P.J. Giver [11] or C. Udri§te [18] agreed that many ap- 
plicative problems coming from Physics [11], [17], Biology [6], [9] or Economics 
[18] can be modelled on 1-jet spaces. 

In such an aplicative-geometrical context, a lot of authors (G.S. Asanov [1], 
D. Saunders [15], A. Vondra [19] and many others) studied the contravariant 
differential geometry of 1-jet spaces. Moreover, proceeding with the geometrical 
studies of G.S. Asanov [1], the second author of this paper has recently elabo- 
rated that so-called the Riemann-Lagrange geometry of 1-jet spaces [5], which 
is a natural extension on 1-jet spaces of the already well known Lagrangian 
geometry of the tangent bundle due to R. Miron and M. Anastasiei [4]. We 
emphasize that the Riemann-Lagrange geometry of the 1-jet spaces allow us to 
regard the solutions of a given GDEs (respectively PDEs) system as horizon- 
tal geodesies [17] (respectively, generalized harmonic maps [7]) in a convenient 
Riemann-Lagrange geometrical structure. In this way, it was given a final so- 
lution for an open problem suggested by H. Poincare [13] (find the geometrical 
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structure which transforms the field lines of a given vector field into geodesies) 
and generalized by C. Udri§te [17] [find the geometric structure which converts 
the solutions of a given first order PDEs system into some harmonic maps). 

In the following, let us present the main geometrical ideas used by C. Udri§te 
in order to solve the open problem of H. Poincare. For more details, please see 
the works [17] and [18]. 

For this purpose, let us consider a Riemannian manifold {M'^,(p^j(x)) and 
let us fix an arbitrary vector field X = on M. Obviously, the vector 

field X produces the first order ODEs system (dynamical system) 

'^^X\x{t)),y i = l,n. (1.1) 

Using the Riemannian metric and its Christoffel symbols 7^-^, and differ- 
entiating the first order ODEs system (jl.ip . after a convenient rearranging of the 
terms involved, C. Udri§te constructed a second order prolongation {single-time 
geometric dynamical system) of the ODEs system (|l.ip . which has the form 

where V is the Levy-Civita connection of the Riemannian manifold (M, (p) and 

is a (1, l)-tensor field which represents the helicity of the vector field X. 

It is easy to see that any solution of class of the first order ODEs system 
(jl.ip is also a solution for the second order ODEs system p.2p . Conversely, this 
statement is not true. 

Remark 1.1 The importance of the second order ODEs system il.S^) comes 
from its equivalency with the Euler- Lagrange equations of that so-called the least 
squares Lagrangian function 

Lis -TM ^M+, 

given by 

Lis{x,y) - \v,,{x) - X^{x)] [y= - X={x)\ . (1.3) 

Note that the field lines of class of the vector field X are the global 
minimum points of the least squares energy action 

¥.is{x{t)) = / LUx''{t),x''{t))dt. 

J a 

As a conclusion, the field lines of class of the vector field X are solutions 
of the Euler- Lagrange equations produced by Lis- Because the Euler- Lagrange 
equations of Lis are exactly the equations U.^) . C. Udri§te claims that the solu- 
tions of class of the first order ODEs system kl.l]) are horizontal geodesies 
on the Riemann- Lagrange manifold [17] 

(RxA/,l + ^,7V(l), =7j,y'=-F;). 
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2 Riemann- Lagrange geometrical background on 
1- jet spaces 

In this Section, we regard the given first order nonhnear ODEs system 
(|l.ip as an ordinary differential system on an 1-jet space J^(T, M"), where 
T C M. Moreover, starting from Udri§te's geometrical ideas, we construct some 
jet Riemann-Lagrange geometrical objects (nonlinear connections, generalized 
Cartan connections, d-torsions, d-curvatures, jet electromagnetic fields and jet 
electromagnetic Yang-Mills energies) which, in our opinion, characterize from a 
geometrical point of view the given nonlinear ODEs system of order one. 

In this direction, let T = [a, 5] C M be a compact interval of the set of real 
numbers and let us consider the jet fibre bundle of order one 

ji(r,M") -^TxR", n>2, 

whose local coordinates i = l,n, transform by the rules 

t = t{t), x' = x'{x'), x\ = ^-^-x' (2.1) 
^ ' ^ ' ^ dxi dt ^ 

Remark 2.1 From a physical point of view, the coordinate t has the physical 
meaning of relativistic time, the coordinate x = {x'^)i=Yn represents the spa- 
tial coordinate and the coordinate y = (2^1 )j=— has the physical meaning 
of direction or relativistic velocity. Thus, the coordinate y is intimately 
connected with the physical concept of anisotropy. 

Let us consider X = ^X^^jj (x'^)^ be an arbitrary d-tensor field on the 1-jet 
space J^(T, R"), whose local components transform by the rules 

_ 9? dt Q-) 

Clearly, the d-tensor field X produces the jet ODEs system of order one {jet 
dynamical system) 

xl=xg(x'=(i)), Vz = M, (2.2) 

where x{t) — (x^(t)) is an unknown curve on R" (i. e., a jet field line of the 
d-tensor field X) and we use the notation 

i not dx^ . - — 

= — — , V 2 = 1, 71. 

^ dt 

Remark 2.2 The main and refined difference between the ODEs systems (OUP 
and 112. which have the same form, consists only in their invariance transfor- 
mation groups. Thus, the ODEs system is invariant under the transfor- 
mation group P = x^{x^) while the ODEs system 112. is invariant under the 
transformation group t — t(t), P = x^{x^) which does not ignore the temporal 
reparametrizations. 
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Now, let us consider the Euclidian structures (T, 1) and (R", Sij), where Sij 
are the Kronecker symbols. Using as a pattern the Udri§te's geometrical ideas, 
we underline that the jet first order ODEs system (|2.2p automatically produces 
the jet least squares Lagrangian function 

JLis : ji(r,M") ^M+, 

expressed by 

n 2 

JL^,(x^a;^)=5][xl-X«(a;)] . (2.3) 

i=l 

Because the global minimum points of the jet least squares energy action 
MiMt)) = J JLis [x\t), — \ dt 

are exactly the solutions of class of the jet first order ODEs system (12. 2[) . 
it follows that the solutions of class of the jet dynamical system (|2.2p verify 
the second order Euler- Lagrange equations produced by JLis {jet geometric 
dinamics), namely 



d[JLu] d fd[JL 



dx^ dt \ dx\ 



As 



= 0, V i = l,n. (2.4) 



In conclusion, because of all the arguments exposed above we believe that 
we may regard the jet least squares Lagrangian function JLia as a natural 
geometrical substitut on the 1-jet space J^(r, R") for the jet first order ODEs 
system ([2^ . 



Remark 2.3 A Riemann- Lagrange geometry on J^(T,M.") produced by the jet 
least squares Lagrangian function JLis, via its second order Euler- Lagrange 
equations (2^, in the sense of nonlinear connection, generalized Cartan con- 
nection, d-torsions, d-curvatures, jet electromagnetic field and jet Yang-Mills 
electromagnetic energy, is now completely done in the works [5], [6] and [7]. 

In this geometric background, we introduce the following concept: 

Definition 2.4 Any geometrical object on the 1-jet space J^(T, R"), which 
is produced by the jet least squares Lagrangian function JLis, via its Euler- 
Lagrange equations i2.4^ , is called geometrical object produced by the jet 
first order ODEs system 1^2. 2\) . 

In this context, we give the following geometrical result (this is proved in the 
works [6] and [9] and, for the multi-time general case, in the paper [7]) which 
characterizes the jet first order ODEs system (12.21) . For all details, the reader 
is invited to consult the book [5]. 
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Theorem 2.5 (i) The canonical nonlinear connection on J'^(T,R") pro- 
duced by the jet first order ODEs system i2. 2\) is 



r = ( 0, N^' 



(i)i 

whose local components N^^^^ are the entries of the matrix 
where 

j(X(i)) = 



is the Jacobian matrix. 

(a) All adapted components of the canonical generalized Cartan con- 
nection CT produced by the jet first order ODEs system 112. 2\) vanish. 



(Hi) The effective adapted components R^^^j^ of the torsion d-tensor T of 
the canonical generalized Cartan connection CT produced by the jet first 
order ODEs system i2.2]) are the entries of the matrices 

Q 

i?(i)fc = ^ [iV(i)] , Vfc = l,n, 

where 

\ ^ / zj — l,n 

(iv) All adapted components of the curvature d-tensor R of the canonical 
generalized Cartan connection CT produced by the jet first order ODEs 
system H2. 2\) vanish. 

(v) The geometric electromagnetic distinguished 2- form produced 
by the jet first order ODEs system ^2.2\) has the expression 



where 



F = fI^^^Sx\ A dx^ , 



Sx\ = dx\ + N^'^^^dx", V i 



l,n, 



and the adapted components F^^^^j o,re the entries of the matrix 



(vi) The adapted components of the geometric electromagnetic d-form 

F produced by the jet first order ODEs system \2. 2\) verify the generalized 
Maxwell equations 
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where ^ ^.-j represents a cyclic sum and 



p(i) ^ 



means the horizontal local covariant derivative produced by the Berwald connec- 
tion BTq on ,/^(T, R"). For more details, please consult [5]. 

(vii) The geometric jet Yang-Mills energy produced by the jet first 
order ODEs system (2. 2\) is given by the formula 

-J n— 1 n ^ 



EYM(x) = - • Trace 



= 1 J=i+1 



In the next Sections, we apply the above jet Riemann-Lagrange geometrical 
results to certain evolution equations from Theoretical Biology that govern two 
of the most actual diseases of our times, namely the spread of cancer cells in 
vivo and the infection by human immunodeficiency virus- type 1 (HIV-1). We 
sincerely hope that our geometrical approach of these evolution equations to 
give useful mathematical informations for biologists. 

Remark 2.6 For more geometrical methods applied to mathematical models 
coming from Theoretical Biology, the reader is invited to consult the book [9J. 



3 Jet Riemann-Lagrange geometry for a cancer 
cell population model in biology 

The mathematical model of cancer cell population, which consists of a two 
dimensional system of ODEs with four parameters, was introduced in 2006 by 
Garner et al. [2]. 

It is well known that cancer cell populations consist of a combination of 
proliferating, quiescent and dead cells that determine tumor growth or can- 
cer spread [3]. Moreover, recent research in cancer progression and treatment 
indicates that many forms of cancer arise from one abnormal cell or a small 
subpopulation of abnormal cells [14]. These cells, which support cancer growth 
and spread are called cancer stem cells (CSCs). Targeting these CSCs is crucial 
because they display many of the same characteristics as healthy stem cells, and 
they have the capacity of initiating new tumors after long periods of remmision. 
The understanding of cancer mechanism could have a significant impact on can- 
cer treatment approaches as it emphasizes the importance of targeting diverse 
cell subpopulations at a specific stage of development. 

The nondimensionalized model introduced by Garner et al. is based on a 
system of Solyanik et al. [16], which starts from the following assumptions: 

1. the cancer cell population consists of proliferating and quiescent (resting) 
cells; 
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2. the cells can lose their ability to divide under certain conditions and then 
transit from the proliferating to the resting state; 

3. resting cells can either return to the proliferating state or die. 

The dynamical system has two state variables, namely P - the number of 
proliferating cells and Q - the number of quiescent cells, and their evolution in 
time is described by the following differential equations {cancer cell population 
flow) 

dP 

--P-P{P + Q) + FiP,Q), 

(3.1) 




-rQ + aPiP + Q) - F{P,Q), 
F(PQ)-J^ r-^ h-^ k-^ 



where 



• a - is a dimensionless constant that measures the relative nutrient uptake 
by resting and proliferating cells; 

• 6 - is the rate of cell division of the proliferating cells; 

• c - depends on the intensity of consumption by proliferating cells and gives 
the magnitude of the rate of cell transition from the proliferating stage to 
the resting stage in per cell per day; 

• d - is the rate of cell death of the resting cells (per day); 

• A - represents the initial rate of increase in the intensity of cell transition 
from the quiescent to proliferating state at small P; 

• A/B - represents the rate of decrease in the intensity of cell transition 
from the quiescent to proliferating state when P becomes larger. 

The Riemann-Lagrange geometrical behavior on the 1-jet space J^{T, M^) of 
the cancer cell population flow is described in the following result: 

Theorem 3.1 (i) The canonical nonlinear connection on J^(r, ]R^) pro- 
duced by the cancer cell population flow i^3.1]) has the local components 



(o,ivj;y , I, J = 1,2, 



where, if 



hQ (1 - kP^) HP 
Fp = ^ and Fn 



(1 + kP^f ^ 1 + kP^ 
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are the first partial derivatives of the function F , then we have 



N, 



(1) 
(1)1 



N, 



(2) 
(1)2 



0, 



^a)2 = -^Si - I [(2« +l)P + aQ-iFp + Fq)] 



{2a + l)P + aQ- 



HQ (1 - /cP2 



HP 



l + kP^ 



(a) All adapted components of the canonical generalized Cartan con- 
nection Cr produced by the cancer cell population flow S3.1\) vanish. 

(Hi) All adapted components of the torsion d-tensor T of the canonical 
generalized Cartan connection CT produced by the cancer cell population 
flow h3.1\) are zero, except 



p(i) 

-"•(1)21 



R 



(1) 

(1)22 



r(2) 
-^(1)12 



i (a - FpQ - Fqq) = ^ (a - Fpq) , 



wher 



2hkPQ (3 - /cP2) h{l- kP^) 

Fpp = . 1 FpQ = 



{i + kp^y 



(l + fcP2) 



dFQQ=Q 



are the second partial derivatives of the function F . 

(iv) All adapted components of the curvature d-tensor R of the canonical 
generalized Cartan connection CT produced by the cancer cell population 
flow ^3.1\l vanish. 

(v) The geometric electromagnetic distinguished 2- form produced 
by the cancer cell population flow i3.1\) has the expression 



Fll^]Sx{Adx\ 



wher 



Sx\ = dx\ + ivf^]fcda;^ V i = 1, 2, 



and the adapted components Fj^i^j, i, j = 1,2, are given by 



^(1)1 ~ -^(2)2 ~ 



f^(l) 

(2)1 

1 

" 2 



-^( i')2 = 2 +l)P + aQ-{Fp + FQ 



{2a+l)P + aQ- 



hQ (1 - fcp2 
(1 + kP^f 



HP 



l + fcP2 



8 



(vi) The biologic geometrical Yang-Mills energy produced by the 
cancer cell population flow I13.1\) is given by the formula 



EYM" 



{2a+l)P + aQ 



hQ (1 - fcp2 

(1 + fcP2)2 



HP 



1 + kP^ 



Proof. We regard the cancer cell population flow (|3.ip as a particular case of 
the jet first order ODEs system (|2.2p on the 1-jet space J^(r, M^), with 



and 



^(1) ("^ I ) 
^(1) I 2; ) 



n = 2, x' = P, x'' = Q 

= x^ -x\x^ +x'^) + F{x\x^) , 
= -rx'^ + ax'^ix^ + x'^) - F{x^,x'^). 



Now, using the Theorem 12.51 and taking into account that we have the Ja- 
cobian matrix 



^(^(1)) = 



1 - 2P - + Fp -P + Fq 
2aP + aQ- Fp -r + aP - Fq 

hQ (1 - kP^) 



\ 



1-2P-Q 
2aP + aQ- 



{1 + kP^f 
hQ (l - kP'^) 



-P 



hP 



[i + kP^Y 



1 

aP- 



kP^ 
hP 



1 



/jP2 j 



we obtain what we were looking for. 



Remark 3.2 (Open problem) The Yang-Mills biologic energetical 
curves of constant level produced by the cancer cell population flow 

^3. which are different by the empty set, are in the plane POQ the curves of 
implicit equations 

hQ (l - kP^) ^ 



Cc 



{2a+l)P + aQ 



(l + fcP2) 



hP 

1 + kP^ 



= 4C, 



where C > 0. For instance, the zero Yang-Mills biologic energetical curve 
produced by the cancer cell population flow ^3.1\) is in the plane POQ 
the graph of a rational function: 



Co-.Q 



P (1 + kP^) [h - (2a +!)(! + kP^ 

fcP2) 



a(l + fcP2) -/i(l 

As a possible opinion, we consider that if the cancer cell population flow does 
not generate any Yang-Mills biologic energies, then it is to be expected that the 
variables P and Q vary along the rational curve Cq . Otherwise, if the cancer 
cell population flow generates an Yang-Mills biologic energy, then it is possible 
that the shapes of the constant Yang-Mills biologic energetical curves Cc to offer 
useful interpretations for biologists. 
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4 The jet Riemann-Lagrange geometry of the 
infection by human immunodeficiency virus 
(HIV-1) evolution model 

It is well known that the major target of HIV infection is a class of lympho- 
cytes, or white blood cells, known as T cells. These cells secrete growth 
and differentiations factors that are required by other cell populations in the 
immune system, and hence these cells are also called "helper T cells". After 
becoming infected, the CDA^ T cells can produce new HIV virus particles (or 
virions) so, in order to model HIV infection it was introduced a population of 
uninfected target cells T, and productively infected cells T* . 

Over the past decade, a number of models have been developed to describe 
the immune system, its interaction with HIV, and the decline in CD4,^ T cells. 
We propose for geometrical investigation a model that incorporates viral pro- 
duction (for more details, please see [8], [12]). This mathematical model of 
infection by HIV-1 relies on the variables T(t) - the population of uninfected 
target cells, T*(t) - the population of productively infected cells, and V{t) - the 
HIV-1 virus, whose evolution in time is given by the HIV-1 flow [12] 

s + (p-d)T-- kVT 

m 

-- kTV - ST* (4.1) 

nST* - cV, 
where 

• s represents the rate at which new T cells are created from sources within 
the body, such as thymus; 

• p is the maximum proliferation rate of T cells; 

• d is the death rate per T cells; 

• S represents the death rate for infected cells T*; 

• TO is the T cells population density at which proliferation shuts off; 

• k is the infection rate; 

• n represents the total number of virions produced by a cell during its 
lifetime; 

• c is the rate of clearance of virions. 

In what follows, we apply our jet Riemann-Lagrange geometrical results to 
the HIV-1 flow (|4.ip regarded on the 1-jet space J^(T, R'^). In this context, we 
obtain: 



dt 

dT* 

IT 

dV_ 

'dt 
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Theorem 4.1 (i) The canonical nonlinear connection on J^(r, M^) pro- 
duced by the HIV-1 flow (^TTp has the local components 



(0,ivg^.), I, J = 1,3, 



where -^(^j^ o,re the entries of the matrix 

, / -kV ~kT 
iV(i) = — \ kV kT-nS 

\ kT -kT + nS 

(ii) All adapted components of the canonical generalized Cartan con- 
nection CT produced by the HIV-1 flow (^TTp vanish. 

(Hi) All adapted components of the torsion d-tensor T of the canonical 
generalized Cartan connection CT produced by the HIV-1 flow (^TTp vanish, 
except the entries of the matrices 




where 



and 

fc/2 
i?(i)3 = ( -fc/2 


(iv) All adapted components of the curvature d-tensor R of the canonical 
generalized Cartan connection CT produced by the HIV-1 flow (^TTp vanish. 

(v) The geometric electromagnetic distinguished 2-form produced 
by the HIV-1 flow i4-i\ ) has the expression 

where 

5x\ = dx\ + ivJ^5feda;^ V i T73, 

and the adapted components F^^^j i,j = 1,3, are the entries of the matrix 

~kV -kT 
kT- nS 

~kT + nS 

(vi) The biologic geometric Yang-Mills energy produced by the HIV- 
1 flow is given by the formula 

^YM."'^-\T,T*,V) = i [k^iV"^ +T'^) + {kT-n5f] . 
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Proof. Consider the HIV-1 flow (|4.ip as a particular case of the jet first order 
ODEs system on the 1-jet space J^iT.W"), with 



n = 3, = T, 



V 



and 



X):Ux\x',x^) = s + [p- d)x^ - ^[x^f - kx^x^ 



"(1) ^ ' ' 

X^'^^{x\x\x^) = kx^x^-Sx^, 



X^^^i^{x ,x ,x-^) = nSx^ — cx"^. 
It follows that we have the Jacobian matrix 



TO 

kV 



-6 kT 

V nS -c J 

In conclusion, using the Theorem 12. 5[ we find the required result. ■ 

Remark 4.2 (Open problem) The Yang-Mills biologic energetical 
surfaces of constant level produced by the HIV-1 flow (^TTp have in 
the system of axis OTT*V the implicit equations 

Ec : k^{V^ + T^) + [kT - n5 f = AC, 

where C > 0. It is obvious that the surfaces Eg" some real or imaginar 
cylinders. Taking into account that the family of conies 

Tc : 2k^T^ + k^V^ - 2kn5T + n'^S^ - 4C = 0, 

which generate the cylinders "Ec, have the matrices 



A = 



2fc2 -knS 
fc2 

~kn6 n'^S^ - 4C 



it follows that their invariants are Aq — k^ (n^(5^ — 8C) , S = 2fc* > and 
I = 3fc^ > 0. As a consequence, we have the following situations: 



1. IfO<C < 



-, then we have the empty set E 



o<c<'- 



2. IfC = 



then the surface T,„_ „2s2 degenerates into the straight line 



T= — 
2k 

V ^0 
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3. If C > , then the surface ^2^2 is a right elliptic cylinder of 

8 > 8 

equation 



^c>^ • ^ — ^2 + -p- = 1' ^* e K, 

where a <b are given by 

0-= TT, > 0= ;= . 

2fc ks/2 
Obviously, it has as axis of symmetry the straight line E„_^j2 . 

There exist possible valuable informations for biologists contained in the 
shapes of the Yang-Mills energetical constant surfaces Sc ? 
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